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16-SEP-1984 7:1 AX/VMS Macro V04-00 Pa 1 
6-SEP-1984 1 99: 93; 8 EMTHRTL. SREIUVXPOWGG.MAR: 1 ” (1) 
TITLE uyxsrouas - G REAL*8 ** G aga. se ower routine 
. IDENT -006/ ; File: OTSPOWGG.MAR EDIT: JCW2006 


'BPARRBARARARARAAAAAAAAALALA AAS ASASEAS ARAL ASE SERRA REALE ALAR EARS SELES AAAS SS 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND sorte 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE pith od NOTICE. Mab SOFTWARE OR ANY OTHER 


* 
® 
& 
® 
* 
& 
* 
o 
0 w 

OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY ®* 
TRANSFERRED. ’ 
e 

be gore? IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
* 

* 

* 

* 

* 

® 

* 

* 


BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 


DIGITAL ASSUMES NO RESPONSIB 


IBILIT 
SOFTWARE ON EQUIPMENT WHICH IS N 


Y FOR THE USE OR RELIABILITY OF ITS 
NOT SUPPLIED BY DIGITAL. 


e@e Ge Ge Se Se Ge Ge Ge Ge Se Se Ge Se Se 
eaeeeeeeeeeaeeeaeeennnnene 


LARRBRAAAAAASSAAAAAREARASALALALALASALASARARARRAAASASESARR ALARA RASA S ALAR AS SD SD 


pf ACILITY: Language support library - user callable 
“ABSTRACT: 
G double base to G double power. 
woes ig eae kh occur if: 


ye base : 
ase and power is 0 or negative. 


VERSION: 2 
HISTORY: 


AUTHOR: 
Bob Hanek 8-Mar-83: Version 2 


MODIFIED BY: 
Jeffrey C. Wiener 16-MAY-83: Version 2-002 


SOOCoCSCOCSCOSCOSCOOSOSOOOCOOOSOSOSOSOSOOSOSOSCOSOSOOOCODOOOOCOO OOOO OOOOOOOoOO 
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; Detailed current edit history 6-SEP-1 
5 -SBTTL HISTORY 


MTHRTL.SRCJUVXPOWGG.MAR; 1 
; Detailed current edit history 


wry 
~ooo 


; Edit history for Version 2 of OTS$POWGG 


5-00) - Implemented new algorithm. RNH 8-Mar-83 
-002 = MicroVAX requires that G_floating routines may not be backed-up 
by awe er a point instructions. To accomplish this the code 
from OTSSPOWDD was modified to create this present version of 
UVXPOWGG. ALL constants had to be scaled accordingly when taken 
from the cAginal OTS$POWGG to presmnents the difference between 
z and 2'=2°6/\n(2). JCW 16-MAY-83. 
3 = Change table INDEX to be local instead of global. LEB 24-May-1983 
4 = Change rererouse of INDEX(Rx) to be INDEXCRx] to avoid Linker 
errors. LEB 25-May-19 
5 = Changed MTH$POWGG entry to OTS$POWGG entry. JCW 26-May-1983 
6 Corrected a bug involving a SYS_F_FLTOVF_F error during a MULG R6, R2. 
Code was added to see if a MTH overflow message or a zero should be 
returned. JCW 19-Jan-1984 


~SBTTL DECLARATIONS 


; INCLUDE FILES: 


EXTERNAL SYMBOLS: 


~DSABL GBL 

-EXTRN MTHSSSIGNAL ; Error signal routine 

~EXTRN = MTHSK_UNDEXP ; Undefined exponentiation code 
-EXTRN MTHSK_FLOUNDMAT ; Floating point underflow code 
-EXTRN MTHSK_FLOOVEMAT ; Floating point overflow code 

> MACROS: 

, SSFDEF ; Define stack frame symbols 

> EQUATED SYMBOLS: 

; base = 4 ; base input formal - by-value 
exp = 12 3; exponent input formal - by-value 
fexp = 8 ; puponent when base is floating 

;_by-value 
ACMASK = “M<IV, R2, R3, R4, RS, Re. R7, RB, RY, R10> 


; OWN STORAGE: None 


PSECT DECLARATIONS: 


—— ~QOODDODO0OO0O0O0O0 00000 0090009 09 09 09 09 09 09 08 SI NIN SINS NIN SII AAA AAAAOAO 
PRO ODN UE WN OOD NAUE WN 9 ODNAULS WN 0 OONOAUE WN O OONOUES WhO 
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-006 DECLARATIONS 6-SEP=1984 :29:1 MTHRTL.SRCJUVXPOWGG.MAR; 1 
0000 113 -PSECT _OTSSCODE PIC,SHR,LONG,EXE,NOWRT 
Wie ; program section for OTS$ code 
138 3; Constants 
8 117; 
a 
00000000 00003CA0 09 ! 9 TWO_M52:.QUAD %*x0000000000003CA0 
000 122 ; 
; The table gives the appropriate entry in the A_TABL 
09 ! The INDEX bl j h i i he A_TABLE 
125 ° 
00 1 INDEX: 
03 Oe Oe 02 01 01 01 69 $60 1 ; -BYTE “X00, “X01, “X01, “X01, “x02, “X02, *x02, *x03 
05 05 05 04 04 04 03 03 001 128 -BYTE “X03, *XO03, *X04, *X04, *X04, “X05, “X05, “x05 
08 07 07 07 07 06 36 88 001 129 -BYTE “X06, “X06, “X06, “X07, “X07, “X07, “x07, *x08 
OA OA OA 09 09 09 08 08 00 130 -BYTE “*xX08, “x08, “x09, “x09, “X09, “XOA, *XOA, “XOA 
OF OF OF OF OF Op OP OD 0080 132 "BYTE “KOD, “OB, “AOD, “KOE, “OES “KOE, “XOEL “KOE 
11 15 1 10 1 OF OF OF $038 138 “BYTE ‘“XOF, “XOF, “XOF, °X10, “X10, *X10, *X10, °X171 
13 12 12 12 12 11:11:11 «40040 134 BYTE SRTh “ite “Nils “mide “Ride “Ride “Ele “eS 
14 14 14 14 14 13 13 13 0048 135 BYTE “X13, “X13, “X13, *X14, *X14, X14, *X14, “X14 
16 16 16 16 15 15 15 1 0050 136 BYTE “X15, “X15, “M15, “X15, “X16, “X16, “X16, *X16 
18 18 18 17 17 171717 #0058 137 -BYTE “X17, “X17, *X17, *X17, *X17, “X18, *X18, *X18 
1A 1A 19 19 19 19 18 18 0060 138 -BYTE “x18, “X18, “X19, *X19, “X19, “X19, “KIA, “XIA 
1B 18 18 18 18 1A 1A 1A 0068 139 BYTE “X1A, “X1A, “X1A, *X1B, “*X1B, “X1B, “X1B, “*X1B 
1D 1D 10 1€ 1€ 1€ 1¢ 1C 0070 140 BYTE PET “Ribas “Ribs “RiGee “Bigs “Rite “Rite “Ete 
1E 1E€ 1€ 1€ 1€ 1D 1D 1D 0078 141 BYTE "RID, “Nits “KIO, “Ris “Riese “Bike “Rieke “HIE 
20 20 20 1F 1F 1F 1F 1F 0080 \*¢ BYTE “Sit, “Ei. “Ht, “HF. *E8, “Eee, “Kee, “hee 
0088 144 
0088 145 ; For k = 0, 2,...., 32, the k-th entry in the A_TABLE is the value 2*(k/32) 
0088 146 : rounded to 113 bits. 
0088 147; 
0088 148 
0088 149 A1_TABLE 
99000009 90006019 08 150 QUAD 240000000 99004010 
8574031 2B00019 444 151 -QUAD “*X8574031 re dt 
B9IOF6CF9 83 84010 009 13 UAD “*X890F6CF9B5584010 
eRe pole 13014011 QOAO 15 QUAD *xX5B851001213014011 
517B3C7D 72 rete OOA8 154 QUAD %*x517B3C7D7 Fete 
igor 4874011 008 155 QUAD “XB9AA 18 4874011 
Ges506/ 8 svt 136 QUAD mtr Leh 
DEEIFSIF 9E9D401 C 1 QUAD “XD F F51F9E9D401 
tae ty pores! oc 138 -QUAD “XB SH A yet 
$768 73A 71A7401 0D 15 QUAD *XA9 8 73A71A ttt 
34 ttt ttre} OD 160 QUAD *X34224C1 pEAGs et 
3 D6 4E£0 rth E 161 QUAD “X89 6 4E084014 
A 3 Aly 14 E 166 QUAD “KX2A g oof A4014 
4F $ 6 4 {Bie F 16 UAD “X4F § 6 4 4 7 15 
54290048 AB07401 F 164 QUAD %*X54290D48AB074015 
5 58034 247E4016 19 165 -QUAD %*X5585B ghee E4016 
0 F 9 16 1 196 QUAD “XxX acpe aes 16 
F 46 EC 1F75401 11 16 -QUAD “X5F af £C1F754017 
0187735EB peta tT 11 168 QUAD “*X ist EBA1144017 
CE13994C 25894018 012 169 QUAD “XCE13994C25894018 
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EAL* G REAL*8 giving 6-SEP-1984 MTHRTL.SRCJUVXPOWGG.MAR; 1 
~SBTTL OTSSPOWGG - G REAL*8 to G REAL*8 giving G REAL*8 result 


<> 


{FUNCTIONAL DESCRIPTION: 
OTSSPOWGG - G REAL*8 result = G REAL*8 base ** G REAL*8 exponent 
The G REAL*8 result is given by: 
base exponent result 


VVvVV A 
ooo Oo O&O 


>0 0.0 
= 0 Undefined Exponentiation 
<0 Undefined Exponentiation 


: any Undefined Exponentiation 
: ¢ gene * Log2(base)) 

: 2 Lexp * Log2(base)) 

: Floating Overflow can occur. 

: Undefined Exponentiation can occur if: 

: ase is 0 and exponent is 0 or negative 
3 2) base is negative 

: CALLING SEQUENCE: 

: power.wg.v = OTSSPOWGG (base.rg.v, exponent.rg.v) 
: INPUT PARAMETERS: 

: base and exponent parameters are call by value 


; IMPLICIT INPUTS: 

; none 

; OUTPUT PARAMETERS: 
; none 


; IMPLICIT OUTPUTS: 
: none 


WWI @ © OO NAULS WIN | 9 ODNAUE WN O OONAUE WH OOONOULS WN Da 


; FUNCTIONAL VALUE: 
: OTSSPOWGG - G REAL*8 base ** G REAL*8 power 


; SIDE EFFECTS: 


SIGNALS 
SIGNALS 


ft ng overt iee 
SIGNALS a 
2) 


: guneer flo Ow if underflow detection is qranies 
P (82 = * UNDEFINED EXPONENTIATION’) if 

~ x” 0 and exponent is 0 or negative 

se is negative 


FUN OOOnNOu 


esSsss 
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WOO OOOO 0000069 09 0909 09 09 09 09 SI INN NNN NN IDS OO DOOD TUT ES ES EE 
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0.0009 00 Cd 0d Co C9 Cd 0d Go G0 CD 0D 
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VXSPOWGG = G REAL*8 * 
2-006 ; 
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Qa 


G 
R 
9 
-ENTRY OTSSPOWGG, ACMASK 3 

¢ ; Move x to RO/RI. If x <0, or x = 0 and y =< 0, return "UNDEFINED 
7 ; EXPONENTIATION’ error condition, otherwise attempt to compute x*ty 
e 

8 


GET_BASE : 
OVG base(AP), RO 3; RO/R1 <== x 


1¢ 3614 BGTR DEF INED : If x > 0 attempt to compute x*tty 
07 (19 BLSS UNDEF INED ; Branch to error code for x < 0 

OC AC 53FD TSTG Ot 3; Check sign of y (Note that x = 0) 
01 3815 BLEQ UNDEF INED ; Branch to error condition if y =< 0 


F A ty eaten, continues here, this implies that x = 0 and y > 0. Return 
3 utty s 
04 RET ; Return 


If yt eat te continues here, this implies that an undefined exponentiation 
was attempted. Signal error and return 


UNDEF INED: 
a. _—e RO 3; RO/R1 <== Reserved operand 


MOVZBL #MTHSK UNDEXP, ~(SP) 
CALLS #1, G*ATHSSSIGNAL 


50 §=—. 8000 He 3¢ 


OO'S8F 9A 


43 Put k 
00000000'GF 01 FB ut error code on stac 


; Convert error number to 32 bit 

3 condition code and signal error. 

; NOTE: Second argument is not re- 

; _ quired since there is no JSB entry. 


RET ; Return 


: If proqeta ing continues here will attempt to coupute xtty as Le a ani 
; We begin by determining k and f such that x = 2*k*f, where 1 =< f < 2. 


DEFINED: 
BICL3 #*XFFFF8O0F, RO, R4 3 R4 <-- 8 deo exponent of x) 
SUBL #*X4010, R4 ; R4 “4e(exponent of x - 1) = 2*4tk 
SUBL R4, RO ; RO <-- f = 2e(fraction field of x) 


; We are now ready to compute log2(x). This computation is based on the 
: following identity: 

: SS eter f-a 

; lLog2(2*k*f) = k + log2(a) + ----- > w--- 2°(2j+1), where z = ----- ° 
3 in(2)/___. 2j+l f+a 

; We begin by determining a as b*i, where b = 2°(1/32) and i is 

3; between 0 and 32 inclusive. Specifically i is chosen by table look-up 
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; in such a fashion as to minimize the magnitude of z. Since log2(a) = i/32 


; we may write 


a poten to INDEX table 


3 log2(x) = k + 1/32 + z*p(z*z). 
EVAL_LOG2: 
ROTL #3, RO R10 
BICL2 #*XFFFFFF80, R10 
Ove INDEXER101, R10 
MOVAW (R10)CR4], =<(SP) 


the final result, 


f and al. 


MOV. A1_TABLECR10], R8 
=(SP) 


XORL3 RY, RI, 
ADDG3 = R8, RO, R2 
R8, RO 
R2, RO, R8 
BICL2 &#*XFFFFO7FF, RO 
; To compute z2 we note 


R4, RO 


We proceed by coseyt ing z = (f-a)/(fta). 
s 


essary to work with a = 
bits respectively of ‘a’. We begin 


R2, R4 
BICL @“XFFEFO7FF, R3, R5 
te R4, R2, R6 


) <== 2°*5(k + 1/32) 


In order to insure the accuracy of 
necessary to compute z to at least 68 bits. Since no 
back up data type is available, we must compute z in fwo parts: z = 21 + 22, 
where 21 is the high 26 bits of z and 22 is the low 5 

to obtain the desired accuracy it is n 
where al and a2 are the high and low 5 
computing (in single precision) 


zi = (f = al)/(f + al) 


Note that f-al can be computed exactly in 53 bits, but f+al any require 54 
bits. The 54 bit can be determined by the exclusive or of the lo 


bits of z. Further, 
al + a2, 


w bits of 


3; RB/RO <== al 
; SP ==> XOR of low bits of al and x 


(This will used to determine the 
bit of f+al.) 


: 54 
; R2/R3 <-- f + al (rounded) 
; RO/R1 <-- f - al (exact) 


; RB/RI <== 21 (26 bits) 


: z= 21 + 22 = (f = al = a2)/(f + al + a2) 
: ==> 22 = (f = al - a2)/(f + al + a2) - 21 


Now let v = f + al + a2 = v1l_+ v2, where v1 and v2 are the high 26 and low 
53 bits of v respectively. Then 

z2 = C(f = al = 2zi"vl) = (a2 + 21*v2)I)/v 
de begin by computing v1 and f - al - 


R4 <== high benguere of f + al 
R4/R5 <== high lo 
R6/R7 <-- w= vl (exact) 
R4/R5 <== zi*vl (exact) 
RO/R1 <-- f = al = zitvl (exact) 


nc 


ngword of f + al (26 bits 
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: $17 : ; Compute v2 and a2 + al*v2 | 
54 FEOO CF4A 7D 0388 418° MOVQ  A2_TABLECR10], R4 ; RG/RS <e= a2 | 
6E FFFEFFFF 8F CA 1 218 BICL #°RFFFEFFFF, (SP) ; Check if w was rounded 
06 13 98 41 BEQL 1$ ; Branch if not rounded 
56 =F C61 gf teed 9A $18 SUBG TWO_M52, R6 ; Correct for rounding error (exact) 
56 4 4OFD AO 419 1$ ADDG R4, R6 3; R6 <-- v2 
56 28 44FD A4 420 MULG RB, R6 3 R6/R7 <H= 21*y2 
56 4 40FD 8 ny, ? 1 ADDG R4, R6 3 R6/R7 <== a2 + 21*v2 
$i ? ; Compute 22 | 
50 56 42FD G AC 425 SUBG R6, Re 3 RO/R1 <-- (feal-z1*v1)-(a2-z1*v2) 
50 52 46FD 0385 426 bived RSL R S RO/RI <-= 22 pit Erg itt | 
0384 428 : The next step is to compute log2(x) accurate to at least 68 bits. This is 
baee : , : accomplished as follows, let | 
0384 431; w= 2egstog2(x) 
0384 +38 3 = (2*5)Ck + 1/32 + z2*p(z*z)) 
0384 433; = 2°5e(k + 1/39) + (2*5)*2"(c0 + c2e2*3 + 1... + €10*2°11) 
0384 434; = “Se(k + 1/32) + z'J] + z'(d2"z'*2 +... + “d10*2 *“10) 
0364 435; = C2*5e(k + 1/32) + 2°) + 2'eQ(z'*2') 
gst tie hell 
0384 438 ; where z' = (2*5*c0)*z and wl and w2 are the high 26 and low = bits of w | 
0384 439; respectively. Note that the choice of ‘a’ used in computing 2. uarantees 
0384 440 ; that.z' overhangs z'*q(z'*z') by at least 15 bits. Hence, if w is computed 
pees ot) ; as wi + we, the necessary 68 bits of accuracy can be obtained. | 
0384 re : We begin by defining | 
0384 445: ¢ = high 53 bits of (2e54c0) 
pyee 446 ; cl = high ¢¢ bits of (2*5*c0) 
beac ret : 1s c2 = low 53 bits of (2*5*c0) 
3 then 
0384 449; z* = (21 + 22)*(cl + 2) 
0384 450; = zitcl + z1*c2 + 22%c. 
gt i 
“yore te i my fot seus sx she 
° : == ¢¢e 
50 FEF2 CE Garp Oset 455 MULG? CRO : RO/RI <= ctze | 
50 58 40FD O3¢7 456 ADDGe R&, RO > RO/R1 <== cez2 + c2e21 
58 50 54 41FD CB O45 ADDG R4, RO, RB 3; RB/RO <== 2° 
0300489 | 
03) £69 > We proceed by letting 
300 462 : wi = high 26 bits of 2*5*(k + i/32) + z1*c1 | 
2 465 ; and | 
4 rr ¢) : w2* = (C2*°5*(k + 1/32) * z1%c1 = wi) + z1*c2) + 22. 
0300 46 : ==> 2*°5e(k + 1/32) + 2" = wl + w2'. : 
0300 468 ; ==> w= t2 Se(k + 1/32) + 2 7 + 2'eq(2"#2") | 
D 469 ; = + we’ + z*aq(z'#2') 
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D 470 ; = wil + w2, 
oe 2 tei dd oad (2"#2") 
; where w2 = w2' + 2'*q(z'*z 
BO ene | 
7E 04 AE are D 475 CVTL 4(SP), =(SP) 3 (SP) <#= 2*5(k+i/32) 
56 4 E 41FD D of8 ADDG (SP), RG, R6 5 RE <m=— 2*5(k 41/32) + z19¢1 = wl 
57 _FFFFO7FF BF CA O3DA 47 BICL2 #*xFFFFOPFF, R7 > R6/R7 <== wl (high 26 bits) 
52 ‘ 43FD E1 478 SUBG ($P)* R6, R2 : R2 <= bits of zi*cl included in wl! 
42FD O366 47 SUBG2 R2, Rb : 
6E 0 4 41FD cf : , ADDG R4, RO, (SP) ; (SP) ==> w2' 
EF 4 ¢ 
EF 4 3 
EF 484 ; Compute w2 
EF 485; 
EF t39 
50 8 8 4S5FD EF 4 MULG3 RB, RB, RO zy RO/R1 <-= z'e2' 
FE9D CF 0 0 55FD O3F4 488 POLYG RO. #LOGLEN-1, LOGTAB : RO/R1 <== g(z'#z") 
50 44FD FB 489 MULG R8, R 3 RO/R1 <W= 2'#Q(z'*2") 
50 6 6E 40FD 0 FF 490 ADDG2 (SP), RO + RO/R1 <== w2 
108 92 
$60 493 ; We now calculate yetogeta? = (ylty2)*(witw2) = ylewl + y2ewl + yew2, where 
rt tee ; yl and y2 are the high and low 26 and 27 bits of y respectively. 
04035 49 
ott 49 
54 OC AC 7D ret: i38 MOVaQ exp(AP), R4 3 R4&/R5 <HW- y 
0407 500 ; Test for the possibility of overflow in the computation of y*yt. 
0407 23) : This wilt occur if the exponent of y plus the exponent of wil is greater 
0407 8 3; than 1023. 
0407 503; 
ot be 
52 54 0B 04 EF 407 05 EXTZV #4, #11, R4, R2 ; biased exp of ¥ 
52 0400 8F A2 O040C 506 SUBW2 #*x400, R2 3 unbiased exp of y 
53 56 OB 04 EF 0411 507 EXTZV #4, #11, R6, R3 : biased exp of w 
53 9300 a - 0416 508 SUBW #*x400, R3 ; unbiased exp of wl 
3 oe 8 509 ADDW R2, ; unbiased exp of wity : 
53 OSFF 8F Bi 041 510 CMPW @*X3FF, R3 ; largest unbiased exp possible is 10235 
03 18 046 11 BGEQ NO_SYS_OVER_FLOW 
0094 31 04 \¢ BRW Y_TIMES_wi_OVER 
4 : 13 NO_SYS_OVER_FLOW: 
5 33 44FD 4 14 MULG2 R4, R ; RO/R1 <-- yove 
5 4 00 ER MOVL R4, R : Rs <-- high longword of y1 
53 55 ware FF 8F cB 42F 216 BICL @*XFFFFO7FF, RS, R3 3 R2/R3 <-- yl 
4 52 42FD 0437 51 SUBG2 R2, R4 > RG/RS <== y2 
52 6 44FD 0438 18 MULG . > R2/R3 <-- ylewi 
54 56 44FD rt! i? MULG R6, R4 > R4&/RS <== y2*wl 
44 2 
7 ¢ ; The next step in computing 2*Cy*log2(x)] is to write ytlog2(x) as 
re 2 : ytlog2(x) = 1 + j/32 + g/32, 
44 6 ; where I is an integer, j is an integer between 0 and 31 inclusive, and 
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rr : ; g9 is a fraction in the interval (-1/2, 1/2) 
44 ; 
56 52 00000000 44494 3§ F G4IFD 044 0 ADDG mths f R2, R6 ; 
56 00000000 00004 F ria 45 1 SUBG #*xX4358, R6 ; 4: <-- 2°5(1 + j/32) 
5 6 42FD 045 § SUBG R6, Re 3 R2/R3 <== those bits gf zityl ngs 
46 ; included in 2°5(1 + j/32) 
28 2 40FD 046 4 ADOG R2, R4 ; 
0 4OFD 0464 5 ADDGe RS, RO > RO/R1 <== g : 
58 8656 waa ¢'¢': = og8 $ BICL #*x00008000, R6, RB 
5 4100 BF Bi 04 CMPW #*X4100, RB 3; Overflow/underflow if exponent field of 
47 8 3 Dok Al te >= 16 
9 19 067 9 BLSS EXCEPTION_1 ; Branch if iI; is tog large 
57 6 4AFD rf A re CVTGL R6, R7 : <-- 2°5#(] + j/32) in integer 
0478 $49 $ 
478 e? ; We can now compute 
0478 44 ; 
Bere et? : xety = 2*Cy#log2(x)J] = 2°CI + j/32 + g/32] 
nore rt! 3 = (2*1)*CA*(B+1)] = 2*I1*CA + A*BJ, where 
0478 549 ; A = 2*(j/32) is obtained by table look-up and B = 2*(g/32) - 1 is obtained 
0478 550 ; by a Min/Max approximation. 
7 
FEGE CF 06 50 55FD 0476 226 POLYG RO, #EXPLEN-1, EXPTAB 3 RO <-- B = 2%*(g/32) - 1. (Chebyshev 
048 554 : polynomial approx. of degree 5) 
52 57  FFFFFFEO 8F CB 04 555 BICL3 #*XFFFFFFEO, R7, R2 3 R2 <-- index into Al_TABLE 
50 FBFB CF42 44FD 068A 556 MULG AI _TABLECR24, RO : RO/R1 <== A®B 
50 FCF9 CF42 40FD 0491 557 ADDG BH » RO : RO/R1 <-- AtB + A2 
50 Fee CF42 4OFD 0498 558 ADDG Al_TABLECR2], R 3: RO/R1 <= 2°C(j+g)/32] = (A*B+A2)+Al 
,.. AA 049F 559 BICW #°R1F, R 3 R7 <== 2°58] 
57 57 FF 8F 9C Q6A 560 ROTL #-1, R7, R7 3 R7 <= 2448] 
50 57 AO OQ4A 561 ADDW2 R7, RO : RO <m= 2*1#2*C(j+g)/32) 
OF 50 81 Q4AA 266 CMPW RO, #*XF 3; test for over/under flow 
08 15 O4AD 56 BLEQ EXCEPTION_2 3; see what exception is if neg or = 0 
04 OQO4AF 564 RETURN: RET 3 otherwise return result in RO 
480 565 
0480 566 ; ee 
vr 200 ; Handlers for software detected over/underflow conditions follow 
0480 569 
48 70 EXCEPTION_1: 
53FD 048 71 TSTG ; if big ARG > 0 gote overflow 
1F 18 048 BGEQ OVER : handler, otherwise go to 
11 0485 57 BRB UNDER 3 underflow handler 
487 are 
487 75 EXCEPTION 2: = 
57 53FD 0487 278 TSTG R7 ; test sign of I; if I <0 
18 18 O4BA 77 BGEQ OVER 3 go to overflow handler 
st a8 
4BC 80 : y*wl would have caused a hardware system floating overflow error. If y<0, 
4BC 1; then we should return a result of 0 since result = 2“(y*(witw2)). Note, 
rts ¢ 3 y can not be zero. 
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48¢ 4 
4BC 5 Y_TIMES_W1_OVER: 
54 53FD 048C § TSTG RS : if y < 0 no overflow is needed 
13. 14 = «04BF BGTR OVER 3 evertiow 
4C1 8 
ee 
rt} 4! ; Underflow; if user has FU set, signal error. Always return 0.0 
4C1 38 
4C1 94 UNDER: 
50 7C¢ rhe 95 CLRQ RO » RO/R1 <-- 0 
0B 04 AD 06 «€E1 04C 36 BBC #6, SFSW_SAVE_PSW(FP), 2s 
rey 9 ; has user enabled Hiya J under flout 
7E OO°8F 9A O04C 98 MOVZBL pee FLOUNDMAT, -(SP) ; Put underflow code 
00000000'GF 01 FB ee 599 CALLS , G*ATHS$SIGNAL ; Convert to MTH$_ FLOUNDMAT. (32- -bit 
4D 600 3; VAX-11 exception code) and 
re 601 3 signal condition 
04 04D o0¢ 2$: RET 3: return 
0404 60 
Rope 604 ; 
404 605 ; Signal floating overflow, return reserved operand, -0.0 
0404 606 ; 
0404 607 
0404 608 OVER: ; else process for overflow 
7E OO°8F 9A 0404 609 MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; Put overflow code on stack 
50 01 OF 79 0408 £610 ASHQ #15, #T, RO ; RO = result = reserved sperand 
040C 611 ; “0.0. RO will be copied t 
04D0C oi¢ i signal mechanism vector (CHF SL _MCH_RO/R1) 
040C 61 3; so can be fixed up by any error 
O4DC §=614 ; handler 
00000000'GF 01 FB QO4DC 615 CALLS #1, G*MTHSSSIGNAL : signal condition 
04 04E3 616 RET 3 return - RO restored from CHFSL_MCH_RO/R1 
04E4 617 
04E4 618 ~ END 
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Al_TABLE 088 R 
A2-TABLE 00196 8 88 
meee = 47F 
BAS = 0 
¢ B8 R 0 
C1 CO R 8 
C2 R 
C 0 FR 0 
DEF INED 000032D R 0 
EVAL_LOG2 38 FR 0 
EXCEPTION_1 0 180 R 0 
EXCEPTION 2 8 487 R 0 
EXP 2 ¢ 
EXPLEN = 00000 
ane se 
INDEX $500008 R 0 
LOGLEN = 00000004 
LOGTAB 00000298 R 02 
MTHS$SIGNAL eeneeeee =X =O 
MTHSK_FLOOVEMAT eeeeeeee xX 00 
MTHSK—FLOUNDMAT teeeeeee =X 00 
MTHSK_UNDE XP ereetaee x 00 
SYS _OVER_FLOW 00000428 R 0 
S$P00 00000308 RG 0 
00000404 R 0 
RETURN 44 4AF R 02 
SF$W_SAVE_PSw s 88 004 
TWwO_AS2 00000000 R 02 
UNDEF INED QOO0031A R 0s 
UNDER 000004C1 R 0 
Y_TIMES_W1_OVER 000004BC R 02 
beeen soe ae eee + 
! Psect synopsis ! 
es 
PSECT name Allocation PSECT No. Attributes 
. . 00000000 0.) 00 ¢ 9. ) NOPIC USR CON ABS LCL oo 8 a NOEXE NORD NOWRT NOVEC BYTE 
SABSS$ 00000000 ( 0.) O1 ¢ -) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_OTSSCODE OOO004GES ( 1252.) O02 ¢ d° ) PIC USR CON REL LEL SHR EXE RD NOWRT NOVEC LONG 
Seqrereninninecapamninmcipaantninanniidbiiitiien ts 
! ; Performance indicators : 
Phase Page faults CPU Time Elapsed Time 
Initialization 2% 68 00:00.1 0:00:00.69 
commend processing 199 0:00: 23% 0:00:07: 
— table sort 0 :00:00.04 0: 00: 4 
Sy Sot table output os a \:29 °: 6:64 
ym able outpu 7:00:00. :00:00. 
Psect synopsis output 3 :00: o 3 : 00: 3 
crene~thterense rie 80; 8: “8 0: §: 5°00 


UVx G 

VAX-11 Macro Run Statistics 

Assembler run totals 400 00:00:04.08 00:00:16.72 

Ing working ssi Limit was 1200 pages. 

1 64 bytes, | pages) of virtual memory were used to buffer the intermediate code. 
re were 


g source Lines were read in Pass 1, producin object records in Pass 2. 
pages of virtual memory were used to define / macros. 


: 
jnqecewequaseeseenauasese 
Macro Library name Macros defined 
“S255SDUAZB:(SYSLIBJSTARLET.MLB;2 0 (itt "'- 
88 GETS were required to define 4 macros. 


There were no errors, warnings or information messages. 


pages of symbol table space al asonee to hold 61 non-local and 2 local symbols. 
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